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ABSTRACT 

We have used the smoothed particle hydrodynamics (SPH) code GADGET-2 to simulate 
the formation of an elliptical galaxy in a group-size cosmological dark matter halo 
with mass Mhaio ^ 3 X lO^2/i-^M0 at z = 0. The use of a stellar population synthesis 
model has allowed us to compute magnitudes, colours and surface brightness profiles. 
We have included a model to follow the growth of a central black hole and we have 
compared the results of simulations with and without feedback from active galactic 
nuclei (AGNs). We have studied the interplay between cold gas accretion and merging 
in the development of galactic morphologies, the link between colour and morphology 
evolution, the effect of AGN feedback on the photometry of early type galaxies, the 
redshift evolution in the properties of quasar hosts, and the impact of AGN winds on 
the chemical enrichment of the intergalactic medium (IGM). 

We have found that the early phases of galaxy formation are driven by the ac- 
cretion of cold filamentary flows, which form a disc galaxy at the centre of the dark 
matter halo. Disc star formation rates in this mode of galaxy growth are about as high 
as the peak star formation rates attained at a later epoch in galaxy mergers. When 
the dark matter halo is sufficiently massive to support the propagation of a stable 
shock, the gas in the filaments is heated to the virial temperature, cold accretion is 
shut down, and the star formation rate begins to decline. 

Mergers transform the spiral galaxy into an elliptical one, but they also reactivate 
star formation by bringing gas into the galaxy. Without a mechanism that removes 
gas from the merger remnants, the galaxy ends up with blue colours, that are atypical 
for its elliptical morphology. We have demonstrated that AGN feedback can solve this 
problem even with a fairly low heating efficiency. Our simulations support a picture 
where AGN feedback is important for quenching star formation in the remnant of 
wet mergers and for moving them to the red sequence. This picture is consistent with 
recent observational results, which suggest that AGN hosts are galaxies in migration 
from the blue cloud to the red sequence on the colour - magnitude diagram. However, 
we have also seen a transition in the properties of AGN hosts from blue and star- 
forming at z ~ 2 to mainly red and dead at 2: ~ 0. Ongoing merging is the primary 
but not the only triggering mechanism for luminous AGN activity. 

Quenching by AGNs is only effective after the cold filaments have dried out, since 
otherwise the galaxy is constantly replenished with gas. AGN feedback also contributes 
to raising the entropy of the hot IGM by removing low entropy tails vulnerable to 
developing cooling flows. We have also demonstrated that AGN winds are potentially 
important for the metal enrichment of the IGM a high redshift. 

Key words: galaxies: formation, evolution, active, cooling flows, intergalactic 
medium 
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1 INTRODUCTION 

The most massive galaxies are elliptical. Their colours and 
spectra are clearly separated f r om those of blue sp irals (e.g. 
iKauffmann et al.l l2003al . l2004l : iBaldrv et al-lliooi ). Most of 



their light comes from old stellar populations ([Thomas et al.l 
Little o r no star formation is present ly going on in 
these galaxies (IKauffmann et al.ll2003al . l2004h . 

Elliptical galaxies also host the most massive black 
holes. The relations that link the black hole mass to 
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bulge pro perties such as luminosity, velocity dispersion 
and mass (|Magorrian et al.lll998l: iFerrarese Merrittll2QQ0l: 



Tremaine et alll2002l : iMarconi k Huntll2003l : iHaring k Rb3 
20041 1 point to an intimate link between the growth of 



black holes and the formati on of spheroids, wh i ch is 
also seen in their c oevolution (|Cattaneo Bernardil [20031 : 
iHopkins et al.|[2007l aV A number of studies have suggested 
that the relation between black hole mass and bulge ve- 
locity dispersion arises because black holes grow until they 



their accretion ( 


Silk & 1 


EleesI 19981: iFabian 1999|: King l2003l: 


Di Matteo et al. 


l2005l: 


Murrav et alJ l2005l: ISpringel et alJ 


2005al: Fabian et al. 2007V This proposal is attractive be- 



cause it also provides a mechanism for terminating star 
formation after a phase of rapid black hole growth (e.g. 
ISpringel et al ] |2005al : iHopkins et alll2007l a.b). 

Previous studies of the role of AGNs in the forma- 
tion of red elliptical galaxies have concentrated on two op- 
posite scenarios. The first one links AGN feedback to the 
growth of super massiv e black holes in galaxy me rgers (e.g. 
Springel et al. 2005a.b: [Hopkins et al.ll2005l . l2007l a.b.c). The 
second one connects the colour bimodalit y to the shutdown 
of cold gas accretion in massive haloes ([Keres et aL |2005[: 
Dekel k BirnboimI [iooi: [Bower et~al[ l2006l: ICattaneo et aU 



20061 . [2007[ )." In the second picture, AGNs play a 'mainte- 



preventing it from cooling teest et al. 2006: Cattaneo et al. 


2006 




Croton et al.[[2006[:[Dunn k Fabian 2006: Fabian et al. 


2006 




Raffertv et al.[[2006[: [Cattaneo k Tevssier|[2007[). Both 



pictures reveal physical aspects of the black hole - galaxy in- 
teraction, but neither of them tells the whole story when it is 
considered in isolation. This is our motivation for studying 
black hole accretion and feedback in cosmological simula- 
tions, where we do not decide for the model at what stages 
of the formation of galaxies black holes should play a role. 

Our first goal is to establish if AGN feedback is re- 
ally necessary to form red elliptical s (see e.g. Birnbo im et 
al. 2007; Khochfar k Ostriker 2007: [Naab et al.[[2007[ ). Our 
second goal is to explore the link between star formation 
and morphology evolution. We want to study how mergers, 
AGNs and the thermal evolution of the IGM contribute to 
the development of the colour - morphology relation. We 
are also interested in the mechanisms that activate black 
hole accretion and in the role of AGN winds in the metal 
enrichment of the IGM. 

We have addressed these questions by simulating the 
galaxy formation process in a dark matter halo that we have 
chosen to be massive enough to form an elliptical galaxy 
(Mhaio — 3 X 10^'^ Mq) without being so massive as to 
make the simulations computationally prohibitive. 

W e have used the GADGET-2 SPH code (|Springel[ 
[2005[ l to integrate the equations of motions for the dark 
matter and the baryons. Sub-grid physics (radiative cool- 
ing, two-phase structure of the ISM, star formation and 
supernovae) are followed w ith the model described by 
[Springel k Hernauist[ (|2003[ ). The model for black hole 
accretion an d feedback is the s a me a s in [Springel et al.[ 
([2005a[ ). The [Bruzual k Charlot[ (|2003[ ) stellar population 
synthesis model has been used for computing magnitudes, 
colours and surface brightness profiles, which we compare 
in simulations with and without AGN feedback. While 
the model of black accretion and feedback that we have 



used has already been employed b y other researchers in 
simulation of isolat e d mergers (e.g. [Springel et al.[ [2005b[ : 
[Hopkins et al.[ [2005[: [Cuadra et al.[ [2006[ ) and cosmological 
simulations (|Li et al. [2007[ v our work differentiates itself 
from previous studies in two respects. The first is that it 
uses high resolution cosmological simulations to investigate 
galaxy formation in a group environment. The second is the 
questions that it addresses, namely the interplay between 
cold flows and mergers in the development of galactic mor- 
phologies, the link between colour and morphology evolu- 
tion, the redshift evolution in the properties of quasar hosts, 
the impact of AGN winds on the chemical enrichment of 
the intergalactic medium, and the effect of AGN feedback 
on the evolution of the photometric properties of early type 
galaxies. The study of the latter has been made possible by 
using a method where not only have we directly analysed the 
simulation outputs, but we have also used them to generate 
virtual astronomical data, which we have then analysed with 
standard observational analysis procedures. 

The plan of the paper is as follows. In Section 2 we de- 
scribe how the simulations are set up both in terms of the 
physical model and of its numerical implementation. Sec- 
tion 3 outlines the methods used to analyse the outputs of 
the simulations. Section 4 and 5 present the formation and 
the evolution of the central galaxy in the simulations respec- 
tively without and with AGN feedback Section 6 discusses 
the effects of black hole heating on the thermal evolution 
of the IGM. Section 7 demonstrates the potential important 
of AGN winds for the early metal enrichment of the IGM. 
Section 8 summarises and discusses the main conclusions of 
this work. 



2 THE SIMULATIONS 

2.1 Initial conditions and numerical parameters 

Using the m ass refinement te chnique described b 



Usmg tne m ass rehnement te cnnique Qescrib e Q by 
Klvpin et aP ([200l[) (see also [Navarro k Bend [l99ll : 
Navarro k Steinmetz[[T997[ ). we have simulated the forma- 
tion of an individual elliptical galaxy within a dark matter 
halo selected for resimulation from a periodic computational 
box of side- length L = 50 h~^Mpc. 

To construct suitable initial conditions, we have first 
created a random realization with N = 2048^ particles 
of a concordance LCDM Universe (^^m = 0.3, Qa = 0.7, 
Qb = 0.045, and h = ifo/lOOkms"^ = 0.7 for the Hubble 
constant). The initial matter power spectrum has been cal- 
culated using numerical results kindly provided by W. Hu, 
and it has been normalized to ag = 0.9. Although as = 0.9 
is not what is currently favored by CMB observations, in the 
both Concordance and WMAP-3 realizations of the Universe 
it is possible to find galaxy with the same properties as we 
have presented. The differences are visible in a statisticall 
sense, but not for a single object. Therefore our choice of a 8 
has no bearing on our conclusions. The initial displacements 
and velocities of the particles have been calculated using all 
waves ranging from the fundamental mode k = 27v/L to the 
Nyquist frequency kny = 27r/L x N^^^ /2. To produce ini- 
tial conditions at lower mass resolution, we have randomly 
selected one in every 8^ (where n=0, 1, 2, 3, 4, 5) neigh- 
bouring particles, assigning to this particle the mass of all 
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the 8^ particles but keeping the position and velocity of the 
selected particle. This procedure allows reducing the mass 
resolution without reducing the small-scale power. By this 
multiple mass technique one can get initial conditions of 
the same realization with different resolutions (e.g. 2048'^, 
1024"^, 512"^ equivalent particle numbers), from which one 
can resimulate exactly the same objects, but with different 
resolutions. 

The initial condition with a resolution of 512^ particles 
has been evolved from z — 50 to 2; = in a simulation with 
dark matter only presented in iMaulbetsch et al.1 (|2007l ). The 
halo finder used to identify dark matter haloes at z = is 
a new MPI+OpenMP version of the Bound Density Max- 
i ma (B DM) algorithm originally introduced bv lKlvpin et al.l 
This algorithm allows to detect isolated or "parent" 
haloes (self-bound structures not contained within larger 
ones) as well as subhaloes (self-bound structures contained 
within larger ones). 

From the halo catalogue identified with the procedure, 
we have selected for resimulation a dark matter halo with 
Mhaio ^ 3.2 X 10^^ Mq. We have considered a sphere of 
radius R — 1.5/i~^Mpc centred on this halo and we have 
used an equivalent resolution of 1024^ particles for the low- 
mass particles inside R. Using several layers of particles with 
progressively increasing masses around the "region of inter- 
est " ensures that the resimulated halo evolves in the proper 
cosmological environment and with the right gravitational 
tidal fields while avoiding large jumps in mass resolution. 
The region of interest is much bigger than the virial radius of 
the selected halo to avoid the penetration of low-resolution 
massive particles into the resimulated object. In addition, 
we have not detected any intrusion of massive particles in 
the region of interest during the whole simulation. 

Finally, we have assumed Qb = 0.045 and we have 
splitted the high-resolution particles into dark matter and 
gas particles. For dark matter particles we have achieved 
a mass resolution of M^^ ^ 8.2 x lO^/i^^M©, while ini- 
tially all baryons are in gas particles with mass M^^^ — 
1.48 X IO^H-^Mq. 

|Pow~ et al.1 (I2OO3I ) give a simple criterion for the gravi- 
tational softening length e necessary in N-body simulations, 
s ^ rvir/V^vir, where rvir and Nvir are the virial ra- 
dius and the number of expected particles within it. They 
find that even a larger softening length keeps the cen- 
tral density profile unaffected. At 2; = 0, our resimulated 
halo has virial radius rvir 300/i~^kpc and virial mass 
Mhaio ^ 3.2 X lO^^/i"^M0. Therefore, for the dark matter 
we have used a gravitational softening length equal to the 
maximum between 0.6 kpc physical, which gives the gravi- 
tational softening length at high and 3/i~^kpc comoving, 
which which gives the gravitational softening length at low 
z. The transition is at z = 4. After z = 4, the force res- 
olution is kept constant in comoving coordinates: 3/i~^kpc 
for the dark matter, 2 h~^kpc for the gas and the stars, and 
l/i~^kpc for the black hole. For the SPH calculations we 
have used 40 db 1 neighbours and we have imposed a mini- 
mum smoothing length equal to 0.75£. 

We have performed two different simulations, with and 
without black hole (see Sections 2.3 and 2.4). The typical 
CPU time for a simulation was ~ 15 CPU days on a SGI 
ALTIX machine with 64 processors. 
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Figure 1. The growth of the stellar mass of the central galaxy 
and of the virial mass of the dark matter halo. The stellar mass 
shown in blue (red) is the one in the simulation without (with) 
AGN feedback. The halo mass (black) is essentially identical in 
the two simulations. The choice of different y-axes for the stellar 
mass (left) and the dark matter mass (right) is such that the 
curves for Mstar and M^aio would overlap if Mstar = /b^halo, 
where /b — 0.15 is the assumed cosmic baryonic fraction. 



2.2 Hydrodynamics, cooling and star formation 

Our initial conditions have been evolved to 2; ^ with an 
updated version of the parallel TreePM-SPH hybrid code 
GADGET-2 (Springel et al. 2001). The code uses an entropy- 
conserving formulation of SPH (Springel & Hernquist 2002) 
which alleviates problems due to numerical over-cooling. 
The code also employs a new algorithm based on the Tree- 
PM method for the N-body calculations which speeds up 
the gravitational force computation significantly compared 
wi th a pure tr ee algorithm. Radiative cooling is computed as 
m iKatz et all (1996) assuming an optically thin primordial 
mix of hydrogen and helium. Photoionisation by an exter- 
nal uniform UV backgro und is computed with the model by 
iHaardt k Madaul (Il996l ). 

The physics of star formation is treated with a subreso- 
lution model in which the gas of the interstellar medium 
(ISM) is described as a two phase medium of hot and 
cold gas (|Yepes et al.] I1997I : ISpringel HernauistI l2003l : 
iMcKee Ostrikerlll977l ). Cold gas clouds are generated due 
to cooling and are the material out of which stars can be 
formed in regions that are sufficiently dense. Supernova feed- 
back heats the hot phase of the ISM and evaporates cold 
clouds, thereby establishing a self-regulation cycle for star 
formation. The heat input due the supernovae also leads 
to a net pressurisation of the ISM, such that its effective 
equation of state becomes stiffer than isothermal. This sta- 
bilises the dense star forming gas in galaxies against fur- 
ther gravitational collapse, and allows converged numeri- 
cal results for star formation even at moderate resolution. 
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Figure 2. The simulation without AGN feedback. Each panel shows a composite UBV image of an 100 X lOO/i^^ kpc^ field. Composite 
UBV images have been generated by combining f/, B and V band .FITS images with the STIFF software by E. Bertin. STIFF makes 
composite colour images by mapping the VBU bands to the RGB colour channels. The snapshots have been selected to illustrate the 
key moments in the evolution of the central galaxy. At z ~ 2.04 — 4.22, the galaxy has the characteristic aspect of a nearly edge-on 
spiral. During the fly-by at z ~ 3.56, a small satellite is tidally disrupted and its debris are accreted by the central galaxy. The galaxy is 
transformed into an elliptical by a major merger at z 1.72 followed by a second major merger at z 2^ 1.28. The visual morphology at 
z ^ 1 is consistent with that of an E+A type object (a red elliptical with a central blue light excess). 



When the SPH density exceeds the star formation threshold 
Pth, the ISM splits into two phases: a cold phase, made of 
star-forming clouds, and a hot phase, made of the tenuous 
gas that fills the space between the clouds. The two phases 
are in pressure equilibrium. They exchange mass and en- 
ergy via cooling, supernova explosions and evaporation of 
clouds by the hot phase. In our simulations, pth corresponds 
to riH — 0.25/i^cm~^. The cold phase makes stars on the 



timescale — t5(p/pth)~"^^^, wh ere tp — 2 . 1 Gmj . These 
values are chosen to reproduce the iKennicuttI (|l998h law. 
Supernova feedback and metal enrichment take place 



^ From a numerical point of view, each gas particle can spawn 
two generations of stars. Therefore, the mass of a stellar particle 
is half of the initial gas-particle mass. 
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Figure 3. The same as Fig. [2] but for the simulation with AGN feedback. The galaxy experienced a major merger associated with 
luminous AGN activity (Lagn ~ ^Edd? ~ —22.6) at z ~ 2.04 followed by another major merger at z ~ 1.72 and a minor merger 
z ~ 1.08 — 1.18. Both produce Seyfert activity at the level of Mb ~ —21.7. The accretion of a tidal feature at 2; ~ 0.63 fuels another flare 
of black hole activity, which is less powerful than the previous ones (M^ ~ ~21) but removes all the cold gas that remains in the galaxy 
causing its transition to the red sequence. A major merger at z ~ 0.075 causes the galaxy to become temporarily bluer in the central 
region. However star formation is quenched by the activation of an AGN (Lagn ~ ^Edd? Mb ~ —23) before any substantial stellar mass 
is formed. At 2; ~ the galaxy has the visual morphology, the photometric properties and the colours of a red elliptical. As the purpose 
of these images is to show the host galaxy evolution, we have assumed that the AGN is obscured (type 2) at all observing times. 



instantaneously whenever stars are formed. For each 1 Mq 
of stars formed, 0.02 Mq of metals are distributed to neigh- 
bouring gas particles. Because we have used a fixed (pri- 
mordial) metallicity for computing cooling, there is no issue 
of using the metallicity of individual particles or smoothing 
with the SPH kernel. Stars have the metallicity of the gas 
out of which they are formed. As in ISpringel HernquistI 



(|2003l ). we adopt a Salpeter (power law) stellar initial mass 
function (IMF) with a slope of —1.3, a minimum stellar 
mass of 0.1 Mq, and a maximum stellar mass of 40 Mq, 
and we assume that a fraction /3sn = 0.1 of the mass that 
forms stars is returned to the hot phase as supernova eject a 
with TsN = 10^ K. The mass evaporated from the cold 
phase of the ISM is larger than this fraction by a factor 



© 2005 RAS, MNRAS 000, ??-?? 



6 A. Khalatyan, A. Cattaneo, M. Schramm^ S. Gottlober, M. Steinmetz, L. Wisotzki 




Figure 4. Density and temperature structure of the IGM in the simulation without AGN feedback (continues on the next page). Each 
image shows a cube of (2h~^ Mpc)"^ centred on the central galaxy. Intensity shows density while hue shows temperature. The white circle 
around the central galaxy shows the virial radius of the dark matter halo. The little subpanel in the lower right corner of each image 
shows an enlargement of the central region. The colour coding is the same as for the large scale image. The histogram on each panel 
shows the density distribution of the gas in the (2h~^ Mpc)^ cube. It spans a logarithmic number density range from 10~^^h^ cm^ to 
lOh^ cm^. The presence of a cold dense phase and a hot dilute phase is clearly visible. The red part of the histogram shows the gas with 
density greater than the star formation threshold tih ~ 0.19/i^cm^. The histogram, therefore, shows the proportion between the masses 
of the hot IGM, the cold IGM and the ISM, which corresponds to the red part of the histogram. At high redshift the gas is cold and 
flows onto the central galaxy in fllamentary streams (z ~ 4.22). At z ^ 2.04 — 2.94, a halo of hot gas has formed inside the virial radius. 
At z ~ 1.08 — 1.42, the cold fllaments fragment and dissolve. By z ^ 0.30 there is little cold gas inside the virial radius. 

oi A = Ao(p/pth)~^^^, where Aq = 10^, so that the hot 2.3 Black hole accretion and feedback 
phase's equilibrium temperature at pth, Tsn/^o, is around 

the maximum of the cooling function. The main difficulty of inserting a supermassive black hole 

in a cosmological hydrodynamic simulation is that we can 
only follow its dynamics adequately when M^h ^ M|^^. 
This occurs because the black hole absorbs the momentum 



© 2005 RAS, MNRAS 000, ??-?? 



Is A GN feedback necessary to form red elliptical galaxies? 7 




Figure 4. Continued. The presence of strong shock due to a sateUite's supersonic motion through the hot IGM and of a tail of ram- 
pressure stripped gas behind the shock is prominent at 2; ~ 0.075. 



of the gas it swallows. At Mbh ^ M|^^, it receives a 'kick' 
whenever it swallows a gas particle. This 'kick' could even 
throw it out of the galaxy. Our solution to this technical 
problem goes as follows. Initially, the black hole has been 
placed at the bottom of the gravitational potential and has 
been made to accrete at a very high rate. As several gas 
particles with momentum oriented at random are accreted 
at the same time, this ensures that the total momentum per 
unit mass accreted is small. We have used this procedure 
to grow a black hole that is essentially at rest with respect 
to the central part of its host galaxy. We have stopped the 
rapid growth when the black hole has reached a mass Mbh ^ 
MS^^ Since in our simulation M^^^ = 1.48 x lO^/i"^M0, 



we have stopped the rapid growth of the seed black hole 
at Mbh = 2 X 10^ h~^MQ. This value is - 10 times larger 
than Mp, but clearly lower than the final black hole mass 
that we expect to find at 2; ^ 0. This prescription forms 
a seed black hole with Mbh = 2 x 10^ Mq already at 
2; :^ 5.4. At that redshift, the host galaxy may not be massive 
enough to host such a massive black hole, although there 
are claims that the black hole - bul^e mass ratio was higher 
at hi^h redshift fe.g lWalter et ahlliooi : iMcLure et alll2006l : 
IPeng et al]|2006l l . We have therefore subsequently frozen the 
growth of the black hole until the host galaxy has reached 
a stellar mass of M^ ~ 1.2 x 10^^h~^ Mq. This happens as 
z 4.2. 
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Figure 5. Density and temperature structure of the IGM in the simulation with AGN feedback (continues on the next page). Colour 
coding and histograms have the same meaning as in Fig. 0] 



Following ISpringel et al.l (|2005af ), the accretion rate posteriori by observing that it produces a reasonable black 

of the black hole has bee n computed with the formula hole mass at the end of the simulation. We have limited 

Mbh = 47ra(GMbh)Vcs"^ (|Bondil Il952l ). The gas density Mhh by requiring that eMbhC^ cannot exceed the Eddington 

p and the speed of sound Cs = ^/(p/lP) have been eval- luminosity LEdd — 1-3 x lO^^(Mbh/lO^M0) ergs~^. 
uated on the SPH kernel of the black hole particle (the Here we have assumed that the power released by the 

40 =b 1 nearest neighbours). The density of the interstel- accretion of gas onto the black hole is eMbhC^ with a stan- 

lar medium smoothed on the scale of the computational dard efficiency of e ^ 0.1, and that a fraction /3 = 0.05 of 

resolution is much lower than the expected high densities this power is converted into heat and distributed to the same 

around the black hole on the scale of the Bondi radius SPH kernel used for computing the accretion rate, 
rsondi = GMbhCs"^ where the Bondi formula applies. We To implement numerically the actual accretion, we fol- 

have compensated for this by deliberately setting a to the low a similar st ochastic approach as i t is applied for regular 

high value a = 300 instead of the value a - 1 (see also star formation (ISpringel et al.ll2005al ). For this purpose, we 

[Hopkins et aL|l2Q06|). The choice of a = 300 is justified a compute for each gas particle j around a BH a probability: 



© 2005 RAS, MNRAS 000, ??-?? 



Is A GN feedback necessary to form red elliptical galaxies? 9 




Figure 5. Continued. 



Pj — ~ — — 5 being swallowed by the black hole. Here, 
Mbh is the black hole accretion rate, At is the time-step, p is 
the gas density estimated at the position of the BH and Wj is 
the kernel weight of the gas particle relative to the BH. We 
then draw random numbers xj uniformly distributed in the 
interval [0, 1] and compare them with the pj. For Xj ^ pj, 
the gas particle is absorbed by the BH, including its mo- 
mentum. On average, this procedure ensures that the BH 
particle accretes the right amount of gas consistent with the 
estimated accretion rate Mbh- 



3 ANALYSIS METHODOLOGY 

Two approaches are possible to extract galaxy masses, star 
formation rates, luminosities and colours from the simula- 
tion outputs. The first is to identify the particle that be- 
long to the central galaxy by running a group-finding algo- 
rithm on the simulation outputs. The second is to convert 
the SPH outputs into virtual observations (.FITS files) by 
using a program that simulates astronomical images. The 
first approach is the one that is commonly used. It has 
the advantage that one has access to the full 6D phase- 
space information in the simul ation outputs. The se cond 
approach has b e en dev eloped bv lContardo et aP (|l998l ) and 
ICattaneo et al.l (|2005al l . It has the advantage that luminosi- 
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Figure 6. Evolution of the star formation rate. The diagram to the left shows the instantaneous star formation rate within the radius 
containing 90% of the galaxy stellar mass at redshift z. The diagram to the right is based on all stars that end up within the 90% mass 
radius at z ~ 0. Therefore, it includes the stars formed in satellite galaxies that merge with central galaxy. The blue (red) lines refer to 
the simulation without (with) AGN feedback. 



ties and colours are determined in a way that is consis- 
tent with standard observational procedures. The two ap- 
proaches complement each other. We have used the first 
one to compute masses and star formation rates, the second 
one to investigate surface brightness profiles and isophotal 
shapes (Sersic indexes, 'boxiness' vs. 'disciness' and colour 
gradients), and both to compute magnitudes and colours. 
The results that we have found with the two methods are 
not quantitatively identical, but the differences are small 
and they do not affect any of the conclusions of this ar- 
ticle. For instance, in the simulation without AGN feed- 
back, the luminosity of the central galaxy estimated by 
the first method is My ^ -21.0, -21.3, -21.8, -20.6 at 
z ^ 0, 0.3, 1, 2, respectively. For comparison, the luminosi- 
ties measured from the .FITS files at the same redshifts are 
Mv ^ -21.1, -21.4, -21.8, -20.4, respectively. The differ- 
ences between the results obtained with the two methods 
are not larger than photometric errors in real data. We now 
expand on the procedures used for these two approaches in 
greater detail. 

3.1 Identifying galaxies in SPH simulations 

To measure galaxy masses, we need a procedure to decide 
which stellar particles belong to a galaxy. For this purpose 
we have used an approximate friends-of-friends group finder 
that links all stellar particles within a linking length equal 
to the force resolution (2/i~^ kpc)0 and have removed all 

The approximate friend-of-friend group finder used 

for this work is publicly available at the URL 

|http: / / www-hpcc.astro.washington.edu / tools / tools.html This 
web site also contains a description of the algorithm. 



particles that are not gravitationally bound to the approx- 
imate friend-of-friend group to which they belong, (notice 
that for a linking length of 2h~^ kpc removing unbound par- 
ticles makes little diff erence to the stellar mass ) . 

We have used the lBruzual CharlotI (|2003[ ) stellar pop- 
ulation synthesis model to compute the [/, B and V band 
magnitudes of the stars inside the galactocentric radius rgo 
containing 90% of the stellar mass Mstar identified with the 
approximate friend-of-friend algorithm. We have made no 
attempt at modelling the extinction of stellar light by dust. 
Let AMstar be the mass of the stars that have galactocen- 
tric distance < rgo and age < At = 10^ yr. We define the 
instantaneous star formation rate as Mstar = AMstar/ At. 

3.2 Observing the simulated Universe 

The approach that we have just described has the advan- 
tage that it uses the full information available in the sim- 
ulations. However, this procedure is quite different from 
the one that an observer would use. Therefore, the values 
that one infers may not be directly comparable to those 
that are measured in observational studies. For this reason, 
we have also devised an alter native procedure. After using 
the iBruzual CharlotI (|2003l l stellar population synthesis 
model to compute [/, B and V magnitudes for each stel- 
lar particle, we have used the software SkyMakeiB (Bertin 
& Fouque 2007) to simulate observations with the Hubble 
Space Telesc ope and convert the SP H outputs into .FITS 
files (see also ICattaneo et al.l (|2005bl )), which has allowed us 
to perform our analysis by using standard software packages 
for astronomical data analysis. 

^ http://terapix.iap.fr/soft 
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Figure 7. The U — V colours of the stellar population inside 
the radius containing 90% of stellar mass. The blue (red) line 
corresponds to the simulation without (with) AGN feedback. The 
vertical grey lines mark the redshifts at which the AGN shines 
with Mb ~ — 21. The colours shown here were computed with 
the procedure described in Section 3.1. Some strong blue colour 
peaks are due to close star-forming companions, which could not 
be disentangled from the central galaxy. 

The analysis of the virtual data has been performed 
on .FITS files that cover a field of 200 x 200/i"^kpc^ 
The field of view is centred on the central galaxy. After 
masking small or distant satellites, we have used GALFIT 
(|Peng et al.ll2002h to fit the surface brightness distribution 
of the central galaxy and its neighbours with a Sersic profile 
E = Eeexp{— 6n[(r/Te)"'^^^ — 1]}, whcrc bn is set by the condi- 
tion that half of the light is contained in the half light radius 
re. By this procedure we have determined the galaxy's mag- 
nitude, half light radius and Sersic index n for each of the 
three bands for which synthetic images have been generated. 

GALFIT has also allowed us to determine the isopho- 
tal shape of our synthetic galaxies. Deviations from purely 
elliptical isophotes are quantified by fitting isophotes with 
the functional form 



Here (xc, yc) define the centre of the ellipses, r the semi- 
major axis and q the axis ratio. For c = 0, Eq. [T] reduces 
to that of an ellipse. The cases of disky and boxy ellipticals 
correspond to c < and c > 0, respectively. 

As our purpose is to concentrate on galaxy evolution, 
even in the simulation with AGN feedback we have assumed 
that the AGN is obscured (type 2) at all times, so that the 
host galaxy photometry is not contaminated by AGN light. 
Therefore, we have avoided the additional uncertainties re- 
lated to the removal of bright nuclear point sources, which 
is needed when studying the host galaxies of type 1 AGNs. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 
z 

Figure 8. Best fit single Sersic index for the major timesteps 
shown in Fig. [2] and [3] The red and blue curves refer to the sim- 
ulations with and without AGN feedback, respectively. Galaxies 
have been classified into discy (round symbols) and boxy (square 
symbol) according to the value of the c GALFIT parameter (Sec- 
tion 3.2) Results for clearly disturbed systems have been marked 
with asterisk and should be taken with caution. The result of the 
simulation with AGN feedback at z = should also be taken with 
caution because we can still see an underlying disc remnant from 
the last major merger. 



4 RESULTS WITHOUT AGN FEEDBACK 

We start by analysing the results of the simulation without 
black hole. The results of the simulation with AGN feedback 
are shown for comparison but they will be discussed in the 
next section. 

The evolution of the luminous component is driven by 
the hierarchical growth of the dark matter halo through 
smooth accretion and merging. The dark matter halo grows 
from Mhaio lO^^h'^MQ at ^ - 5 to Mhaio 3 x 
10^^h~^ Mq at z ^ (Fig. [1]). Two major mergers with 
other haloes are recognisable, the first at z ^ 2.7 — 2.9, the 
second at z ^ 0.7 — 0.9. Qualitatively, the growth of the cen- 
tral galaxy follows that of the dark matter halo. E.g., the 
halo merger at 2.7 — 2.9 is followed by a major merger 
between galaxies at z ^ 2. This merger is clearly visible in 
Fig. [Hand Fig. El 

Fig. [2] and Fig. [3] summarise the key moments in the 
evolution of the central galaxy from 2;~4.2toz^0in the 
cases without and with AGN feedback, respectively. 

Our simulations support the notion that the forma- 
tion of galaxies at the centres of dark m atter haloes is due 
to accretion of co ld filamentary flows (|Keres et al.l l2005l : 
iDekel BirnboimllioOG ). At z = 4.22, the accretion of cold 
streams onto the central galaxy is clearly visible in the en- 
largements of the central region (Figs |4] and [5]). In our sim- 
ulations, the accretion of cold flows form galaxies with disc- 
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Figure 9. U — V galaxy colour profiles in the simulations without (left panel) and with (right panel) AGN feedback for certain key 
moments in their evolution. The profiles without AGN always show the same behaviour with a clear drop toward the centre. In the 
case with AGN feedback, we clearly see standard colour gradients (the colour becomes redder moving inward) after a period of passive 
evolution and reversed colour gradients when star formation is reactivated by mergers. 



like morphologies. An edge-on disc is clearly recognisable in 
Fig. [2] at z 2.23 -4.22. 

The galaxy mass grows proportionally to the halo mass 
until z ^ 4.3. At this redshift, the stellar mass of the central 
galaxy is Mstar ^ 0.7/bMhaio , where /b — 0.15 is the cosmic 
baryonic fraction (Fig. [T]). The situation changes at z ^ 4. 
Due to the higher halo mass, a large fraction of the gas that 
flows into the virial radius is heated by shocks before it can 
reach the central galaxy. Fig.|4]makes this point very clearly. 
From z ^ 2.94, the ^ 10^ K hot spherical halo (in green) fills 
an increasingly large fraction of the volume inside the virial 
radius (white circle). The cold gas that manages to stream 
deep into the halo fragments into clumps before it can be 
accreted onto the central galaxy (see the enlargements of 
the central region). Therefore, already at z ^ 4, the supply 
of cold gas to the central galaxy has started to dry out. 
This is the reason why the growth of the galaxy stellar mass 
between z ^ 4 and z ^ 2 is very limited (Fig. [1]). As soon 
as the galaxy stops to accrete gas, the star formation rate 
begins to decrease (Fig. [6]) and the colour starts to become 
redder (Fig. [T]). This can also be appreciated by comparing 
the colours at z ^ 4.22 and z 2.23 in Fig. [2] Without 
refuelling by galaxy mergers, the galaxy is already on its 
way to becoming an early- type spiral. 

However, while the star formation rate in the central 
galaxy is already steeply declining soon after z ^ 4, the 
total star formation rate within the halo is only exhibiting 
the same steep decline after z ^ 2. Fig. [6] shows this point 
by comparing the central galaxy's star formation rate (left) 
with the cumulative star formation rate of the central galaxy 
and the satellite galaxies that will merge with the central 
galaxy before the end of the simulation at ^ ^ (right). 
The different behaviour in the two panels of Fig. [6] is due 
to the collapse of gas in substructures along the filaments. 
They contribute to the growth of the central galaxy in two 
ways: directly, by depositing gas onto it, and indirectly, by 



forming satellites, which then fall into the central galaxy. 
The second mode becomes the predominant mode of galaxy 
growth after the formation of a hot halo. Comparing the two 
panels in Fig.[6]shows that that the second mode contributes 
to a substantial fraction of the final galaxy mass. 

Another process that becomes important and affects the 
evolution of satellite galaxies after the formation of a halo of 
hot gas is ram pressure stripping. The effects of ram pressure 
stripping are clearly recognisable at z ~ 0.075 and z ~ in 
Fig. m The green 'tongues' to the left of the hot central 
region are trail of gas stripped from galaxies that have en- 
tered the virial radius and moved through the hot gas. In 
front of these 'tongues' strong shocks (red) due the galaxies' 
supersonic motions are clearly visible. 

The quiescent merging history at high redshift is the 
reason why the galaxy has been able to maintain a spiral 
morphology until z 1.96. The interaction at z ^ 3.56 leads 
to the accretion of tidally disrupted material. It causes a 
temporary increase in the star formation rate (Fig. [6]) but it 
has no effect on the galaxy mass. By z ~ 2.04, the galaxy has 
become part of small group. An approximately equal mass 
merger with another disc galaxy at z ^ 1.72 transforms our 
spiral galaxy into an object whose visual morphology is more 
consistent with an elliptical galaxy (Fig.[T]and Fig. [2]). How- 
ever, a quantitative analysis of its photometry shows that its 
Sersic index is more consistent with a spiral galaxy both be- 
fore (n ~ 0.7) and after (n ^ 1.3) the merger. This merger is 
followed by another major merger, with a 1 : 2 mass ratio, at 
z ~ 1.28. The latter occurs after two satellite galaxies have 
first merged into one. After this second merger, the Sersic 
index rises to n 1.6, but is still outside the range that we 
expect for an elliptical galaxy. These two major mergers are 
clearly visible as star formation rate peaks in Fig. [6l The 
peak at 2; 0.8 is due to the accretion of gas stripped from 
tidally disrupted companions. Only after the star formation 
rate starts to decline {z ^ 0.6) does the Sersic index reaches 
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values more typical of an elliptical galaxy (e.g. n ^ 2.7 at 
z ^ 0.3). The accretion of tidal features, one of which is also 
visible in Fig. [6] at 2; :^ 0.075, temporarily reactivates star 
formation by bringing fresh gas into the galaxy. This effect 
is clearly visible as a central light excess in an otherwise red 
elliptical galaxy (Fig. [2]). This picture is also supported by 
the colour gradients in Fig. [9l where we see a steep drop of 
the U — V colour toward the galactic centre Consequently, 
the galaxy colour is not as red as it is expected for a red 
sequence galaxy (Fig. [71 and [TOl) . This find ing i s consistent 
with p revious studies by iMeza et al.l (|2003l l and iNaab et"al] 
([2003), who also reported values oi U — B that were not 
increasing toward the centre. Elliptical galaxies with blue 
central regions are observed (EH- A galaxies), but they are 
not typical of the elliptical population. While an individual 
object cannot be used to draw statistical conclusions, the 
simulations carried out until now support the conclusion 
that a mechanism that removes the cold gas accumulated 
in galaxy mergers is necessary to form truly red elliptical 
galaxies. Elliptical galaxies formed in simulations that do 
not contain such a mechanism have colours in the 'green 
valley' that separates the blue and the red population on 
the colour - magnitude diagram (Fig. [TO)) . 



5 THE SIMULATION WITH AGN FEEDBACK 
5.1 Dynamical evolution 

The dynamical evolution in the simulation with AGN feed- 
back is very similar to that in the simulation without black 
hole until AGN feedback is activated at 2; 4.2. Small dif- 
ferences in the dynamics before 2; :^ 4.2 are due not only to 
the fact that a small mass of gas has been used to form a 
supermassive black hole (Section 2) but also to differences 
in the timesteps of the two simulations. On large scales, the 
dark matter distribution is very similar with and without 
AGN feedback even after z ~ 4.2. Inside the virial radius, we 
continue observe a similar mass distribution of dark matter 
substructures, but substructures are not found at the same 
positions because even small differences become very large 
in a few dynamical times (the evolution is highly non-linear 
and phase mixing occurs). 

The reason why we have not turned AGN feedback 
on from the beginning is computational rather than phys- 
ical (Section 2). Al z ^ 4.2, the galaxy stellar mass is 
lO^°/i"^M0 and the black hole mass is Mbh — 
2 X 10^ Mq. The presence of a large mass of cold gas 
that has accumulated unhindered inside the central galaxy 
means that the accretion rate ai z ^ 4.2 is very large. The 
black hole swallows 5 x 10^ Mq of gas in a very short 
span of time (Fig. [TT]). The energy injection that accompa- 
nies the fast growth of the black hole at z ^ 3.9 — 4.2 causes 
an explosion, which reheats or blows away nearly all the cold 
gas that is present in the galaxy at that time. Star forma- 
tion is quenched until z ^ 3. The timing and the intensity 
of this accretion episode are artificial as they depend on the 
arbitrary choice of the redshift at which AGN feedback is 
activated. However, it is energy and not power that mat- 
ters when we are concerned with the effects on cosmological 
timescales as opposed to the effects on the timescale of the 
active phase, lO'^ - 10^ yr. The total mass of 3 x lO^M© 
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Figure 10. The U — V colour and the V-band magnitude of the 
central galaxy at 2; ~ in the simulation without (blue symbol) 
and the simulation with (red symbol) AGN feedback. The results 
of the simulations have been superimposed on a synthesised U — V 
vs. My colour - magnitude diagram for a volume- weighted sample 
of 14.5 ^ r ^ 16.5 Sl oan Distal Sky S urvey galaxies, which we 
have reproduced from lBell et al.l (|2004l ). 



that our black hole has accreted by 2; 2.6 — 3 fFig. llip and 
thus the total energy that it has deposited into the IGM by 
that redshift are not unreasonable. We also note that the 
stellar mass (Fig. [T]) and the star formation rate (Fig. [6| 
at z ^ 2.6 are similar in the simulations with and without 
AGN feedback. Therefore, we conclude that, although the 



evolution at 3 , 



: 4 is affected by the arbitrary way with 



which we have switched on AGN feedback, the memory of 
the explosive growth ai z ^ 4.2 is lost hy z ^ 2.6 — 3 (Fig. [6]). 
For this reason, we concentrate on the evolution after ^ ^ 3. 

The reactivation of star formation at 2; 3 is due to 
a gas-rich minor merger (Fig. (3] 2; ~ 2.94). Despite the im- 
pression of an elliptical morphology, the remnant at 2; ^ 2.23 
is fitted by a Sersic profile with n c^i 0.7. This merger is 
followed by two major mergers at z ^ 2.04 and z :^ 1.72 
(Fig. (3]). The first of the two is a nearly equal mass merger. 
Their combined mass growth is similar to the one produced 
by the single major merger observed at 2; ~ 1.72 in the 
simulations without AGN feedback. It is interesting to see 
how the different kinematics of the mergers at 2; :^ 2.04 and 
z ^ 1.72 affect black hole growth and star formation. The 
first merger is nearly head-on. Therefore, most of the gas 
in the companion is brought straight into the centre. The 
result is a fiare of AGN activity with Mb ^ —22.6, where 
the black hole reaches the Eddington limit. We also observe 
a period of enhanced star formation clearly visible as a peak 
in U — V colour in FigU] In the second merger, the two 
galaxies merge with a spiral motion (Fig. [3]). Their tidal in- 
teraction induces a burst of star formation before their cores 
have time to merge. It is at that point that the growth of 
the black hole is reactivated and quenches star formation 
(Fig.[6|). In this second case, the peak star formation rate is 
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Figure 11. The growth of the black hole in the simulation with 
AGN feedback. Besides the initial growth, four main growth 
episodes are identifiable. They coincide with the major merger 
at z ~ 2.2, the borderline merger at z ~ 1.1, the accretion of a 
tidal feature at z ~ 0.6 and the major merger at 2; ~ 0.1. Only the 
first and the last correspond to accretion rates approaching the 
Eddington limit (see Fig. [6] and the caption to Fig|3}. Notice that 
this diagram uses a linear rather than logarithmic mass scale. 



higher but the peak AGN luminosity is lower [Mb ^ —21.7) 
because the AGN only switches on in the final stage of the 
merger, when there is little gas left. After the z :^ 1.72 
merger, at 2; 1.42, the galaxy has the appearance of a red- 
dish mildly boxy elliptical galaxy with Sersic index n ^ 2.5. 
The 'noisy' star formation rate history at 2; ^ 0.6— 1.3 is the 
result of a minor merger at 2; 1.08 — 1.18, which activates 
an AGN with Mb ^ —21.7, and of the accretion of tidally 
disrupted clouds from a number of close small companions. 
The gas brought in by these events is the reason why the 
isophotes return to be temporarily discy. The accretion of a 
tail of tidal debris (the blue half ring in Fig. [3]) triggers an 
AGN with Mb ^ —21 at z ^ 0.63. Star formation is, there- 
fore, first reactivated, and then quenched. From this point 
onwards, the specific star formation rate is the same range 
as the values inferred for red sequence galaxies in the Sloan 
Digital Sky Survey (Kauffmann et al. 2004). 

By z ^ 0.1, when a new merger is about to occur, the 
galaxy has evolved into a boxy elliptical with Sersic index 
n ^ 6.8. The increase in the value of the Sersic index and 
the transition from discy to boxy isophotes occur gradually 
and do not appear as immediate results of galaxy merging, 
suggesting that interactions with small companions and in- 
ternal dynamical evolution have played a role in this trans- 
formation. However, another important element is the fad- 
ing of discs of young stars formed by dissipational infall in 
galaxy mergers. These discs contribute to a small fraction 
of the mass of an elliptical galaxy but to a less small frac- 
tion of its light. They therefore contribute to the discy mor- 
phology of the remnants of dissipational mergers. However, 



when the newly formed stars age, their luminosity contrast 
with respect to the preexisting stellar population weakens. 
Therefore, the galaxy evolves into an elliptical with boxy 
isophotes. We have verified this interpretation by consider- 
ing the central galaxy at z ~ 0.91 and aging it passively to 
z ^ without changing anything of its mass distribution 
and kinematics. The galaxy passed from being a discy ellip- 
tical with n 3.1 to being a boxy elliptical with n ^ 4.3, 
although passive fading alone does not bring the Sersic index 
up to n ^ 6. 

A final major merger with a disc galaxy at z 0.075 
is not able to reactivate star formation effectively because 
its gas content is immediately blown away by a powerful 
AGN outbursts with Mb —23. Late mergers are less ef- 
fective at bringing gas into the centre not only because gas 
is consumed by star formation but also due to ram pres- 
sure stripping of the ISM by the hot gas (Section 4). The 
z ^ remnant of this final merger has a Sersic index of 
n ~ 4.0 and slightly discy isophotes related to the presence 
of an underlying disc from the last major merger. Due to 
quenching of star formation by AGN feedback, the galaxy is 
now clearly on the red sequence of the colour - magnitude 
diagram fFig. [TO)) . 



6 THE ENTROPY OF THE INTERGALACTIC 
MEDIUM 

In our simulation, the primary effect of AGN feedback is 
to quench star formation by expelling the galaxy's ISM. 
However, the development of a cooling flow could reactivate 
star formation and move elliptical galaxies back to the blue 
cloud. To address the 'maintenance' problem in the context 
of our simulations, we first need to determine the cooling 
time of the hot gas. 

For this purpose, it is useful to introduce the entropy 
constant S = kTn~^^^ because it allows to writ e the cooling 
time as the product of two terms ([Scannapieco Oh..20o3 ) : 



fn/cT 



Me 



2A(T,Z) VlOkeVcm 



(lOkeV cm^)2. 



(2) 



My (A:T)2A(T,Z) 

The first term is a measure of the entropy and, therefore, 
it does not change when the gas is compressed or expands 
adiabatically. The second term depends on the temperature 
and the metallicity only. For a fully ionised gas, /i ^ 0.59 and 
/ie — 1-1 • In the range of temperatures that we encounter in 
our simulations (T < 10^ K), the factor on the second line 
of Eq. (1) has an absolute minimum. Its value is ^ 2 Gyr for 
the primordial cooling function used in our simulations, and 
^ 200 Myr for Z — O.SZq. Therefore, the entropy constant S 
determines the minimum cooling time for a given metallicity. 
We use the notion of minimum cooling time to interpret the 
evolution of the IGM entropy distribution, where by IGM we 
mean the gas below the density threshold for star formation, 
which separates the IGM from the ISM. 

The IGM entropy distribution is strongly bimodal both 
without and with AGN feedback (Fig. [12]). The low en- 
tropy peak corresponds to the cold gas in the filaments. 
The high entropy peak corresponds to dilute hot gas at 



© 2005 RAS, MNRAS 000, ??-?? 



Is A GN feedback necessary to form red elliptical galaxies? 15 





Log S (keV cm') Log S (keV cm') 

Figure 12. Entropy distribution of the IGM inside the virial radius in the simulation without (blue) and the simulation with (red) AGN 
feedback. The ISM (the gas with p > pth: where pth is density threshold for star formation) has not been shown on these diagrams. The 
vertical dotted lines mark the values of the entropy constant S that correspond to a minimum cooling time of 0.1 Gyr, 1 Gyr and lOGyr 
for primordial metal abundances. 
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the virial temperature (Tvir ~ 2 — 3 x 10^ K at 2; 0) 
and it is where we expect it from the simple calculation 
S = A;Tvir(^^bpvir//imp)"^/^ oc MX^p'.y^. Gas that flows 
into the virial radius at lower z is shocked to higher entropy 
because its density is lower and the potential well of the 
dark matter is deeper. If we compare Fig. [12] with the den- 
sity histograms in Figs. HHSl the gas in the low entropy peak 
corresponds the dense (cold) IGM, while the gas in the high 
entropy peak corresponds to the low density phase of the 
IGM. The ISM (the red area under the density histograms 
in Figs. SHSI) is not considered in Fig. [12] 

As we move from high to low redshift, the cold IGM 
progressively disappears. It contains a small fraction of the 
mass at z ^ 0.3 — 1. By z ^ 0.1 it has almost completely 
vanished. Most of the gas is in the hot IGM, and the cold 
gas that remains is in the ISM (Figs.|4]-[5]). In addition, more 
and more baryons pass the entropy threshold above which 
their cooling time is longer than the age of the Universe. At 



; 0.3 most of the IGM has tc 



10 Gyr. 



Fig. 1 121 shows that AGN feedback does not make a huge 
difference. We now explain how we can understand this re- 
sult by considering the total energy output of supernova 
and AGN feedback integrated over the lifetime of the Uni- 
verse. Our model of supernova feedback assumes that super- 
novae return 10% of the mass that makes stars to the ISM. 
This mass is returned with a temperature of Tsn = 10^ K. 
The energy efficiency of supernova heating is, therefore, 
csn' = f /cTsN(O.lM0//imp)/(M0C^) - 2 x 10"^ The en- 
ergy efficiency of AGN heating is much higher (cbi?* — /3e ^ 
5 X 10~^) but its effectiveness is suppressed by a factor equal 
to the black hole-to- stellar mass ratio at 2; = 0, Mhh/Mstar 
5 X 10~^. When we multiply /3e by this factor, we find an ef- 
fective efficiency of 2.5 x 10~^. Therefore, in our simulation, 
black hole and supernova heating are comparable. This is 
partly because our black hole is not very massive. Its mass 
at 2; = is only Mbh ^ 5.2 x lO'^/i~^M0. Typical black hole 
masses in elliptical galaxies with Mstar 10^^h~^ Mq are 
3 — 4 times larger. However, there is also a compensation 
mechanism. Without AGN heating, there is more star for- 
mation and, therefore, more supernova heating. However, 
strong supernova feedback is only possible in blue (star- 
forming) galaxies. Therefore, the feedback from type II su- 
pernovae included in our simulations cannot be used to pre- 
vent th^_£ea£tivati^^ formation in elliptical galaxies 
(see ICiotti Ostrikerl l2007l for a discussion of the role of 
type I supernovae in this context). 

A closer look at the red and the blue curves in Fig. [T2l 
shows two important differences. First, at 2; ^ 0.075 and z ^ 
0, there is a small but clearly visible tail with t™^^ < 1 Gyr 
under the blue curve. There is no such tail under the red 
curve, which is slightly shifted and skewed towards higher 
value of 5*. In other words, the mean entropy may be similar, 
but the chances that some hot gas will manage to cool in the 
next billion years are much smaller in the simulation with 
AGN feedback. Second, the total mass under the red curve 
is lower than the total mass under the blue curve. This is 
due to the fact that a quarter of the baryons in the hot IGM 
have escaped the virial radius in the case of the simulation 
with AGN feedback. 

In our simulations the reactivation of star formation 
and black hole accretion is mostly due to mergers. Cool- 
ing flows onto the central galaxy do not play a major role 



in replenishing th e central galaxy with cold gas. Work by 
iNaab et ahl (|2007l l suggests that this would remain true even 
if both AGN and supernova fee dback were t u rned off. How- 
ever, this occurs because both Naa b et"al] (|2007l ) and we 
have computed cooUng for primordial metal abundances. 
For Z — O.3Z0, the minimum cooling time is over ten 
times shorter. In fact, there is evidence that feedback from 



solving the cooling 


flow problem (Fabian et al. 20021, 


20031; 


Forman et all 12005 


: Voit & Donahue. 2005l: 


Best et al.l 


200d; 


Dunn & FabianI I2OO6I: iFabian et all I2OO6 


; iRaffertv et al. 


20061), even if these observations are urimarilv concerned 



with galaxy clusters rather than with smaller systems, such 
as the one that we have simulated. 



7 METAL ENRICHMENT 

AGN winds contribute to the metal enrichment of the IGM 
by lifting highly enriched gas from the galactic centre. Their 
potential importance is demonstrated by the diagram in 
Fig. I13[ which shows the mean IGM metallicity as a func- 
tion of its overdensity with respect to the mean baryonic 
density of the Universe. Without AGN feedback, metals are 
conflned to the vicinities of the galaxies in which they are 
produced (blue curve). In contrast, AGN feedback leads to 
a strong increase of the global metallicity of the IGM, which 
is easily overpredicted already at 2; ~ 3 (red curve) . 

The model of the ISM that we have used for this work 
does not allow the effective generation of supernova winds 
because it is based on the assumption that the cold and the 
hot ph ase are in pressure equilibrium. ISpringel HernquistI 
(I2OO3I ) have developed a modified version that does it. They 
have demonstrated that supernova driven winds can also be 
important for the metal enrichment of the IGM, although 
they are not sufficient to account for it quantitatively. 

Fig. [13] shows that the AGN feedback model and the 
efficiency assumed in this study over-enrich the IGM. How- 
ever, what this work demonstrates is that even moderate 
AGN winds could play an important role in the chemical 
evolution of the IGM at high redshift. 



8 DISCUSSION AND CONCLUSION 

We have used cosmological hydrodynamic simulations with 
the SPH code GADGET-2 to study the formation of the 
central galaxy in a halo with a final mass of Mhaio — 
3 X 10^^ Mq at 2; = 0. The simulations include sev- 
eral baryonic processes, e.g. the presence of a photoionising 
UV background, radiative cooling, the formation of a two- 
phase ISM, star formation, supernova feedback, the evap- 
oration of cold clouds and the chemical enrichment of the 
ISM. The model used to des cribe these processes is that by 
ISpringel HernquistI (|2003l V We have also run simulations 
where we inclu de the growth of a central black hole (com- 
puted with the iBondii n 



I952I model) and its effects on the 



ISM/IGM (5% of the power released by black hole accre- 
tion is instantaneously thermalised in the surrounding gas) . 
This mo del for black hole accr etion and feedback was devel- 
oped bv lSpringel et al and has been used in several 
publications (e.g. Springel et al.|[20'05bl : iHopkins et al]|2005l : 
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Figure 13. Mean gas metallicity at z ~ 3 as a function of its 
overdensity with respect to the mean baryonic density of the 
Universe. The blue and red curves correspond to the simulations 
respectively without and with AGN feedback. The grey shaded 
area shows iSchave et al 1 (|2Q03l Vs determination of the metallicity 
of the IGM at 2; ~ 3 as traced by CIV. Although the data them- 
selveses are considerably uncertain, we clearly see that the model 
without AGN feedback seriously underestimates the metallicity 
of the IGM at baryonic overdensities ^ 10, while the model with 
AGN feedback overestimates it. 

ICuadra et aLll2006l : iLi et al.ll2007l ). The main difference be- 
tween this work and previous studies is not in the physics, 
but rather in the system that we have chosen to simulate 
(the central galaxy of a group instead of an isolated merger 
or a cluster) and in the questions that we have addressed. 
This paper is focussed on the the interplay between gas ac- 
cretion and merging in the development of galactic mor- 
phologies, and on the role of AGN feedback in the origin of 
the colour - morphology relation. It is also the first work of 
this kind to discuss the effect of AGN feedback on the pho- 
tometry of early type galaxies, the redshift evolution in the 
properties of quasar hosts, and the impact of AGN winds on 
the chemical enrichment of the intergalactic medium (IGM) . 

Our work supports the notion that the early phases 
of the formation of galaxies are driven by the accretion 
of cold filamentary flows. Our simulations show that this 
mode of galaxy formation only forms discs. Elliptical mor- 
phologies only appear after galaxy mergers have occurred. 
There is a clear parallel between the growth of the dark mat- 
ter halo and the growth of the central galaxy via merging 
(Fig. [1]) . These results were largely expected but they are 
important because they provide a powerful confirmation of 
the same assumptions th at are used in sem i analytic models 
of g alaxy formation fe.g. ' Bower et al.ll2006l : ICattaneo et al.l 
[2OO6: Croton et al. 2006). 

Disc star formation rates at 2; 4 are about as high 
as the peak star formation rates in major mergers at later 
redshifts (in our simulations mergers start being important 
at z ^ 2). That is because the merging galaxies have a gas- 
to-stellar mass ratio of ^ 0.1 by the time the first major 
merger occurs. However, already ai z ^ 4 the filaments 
start to dry out, they fragment and they are no longer con- 
nected to the disc of the central galaxy. Therefore, the cen- 



tral galaxy's star formation rate begins to decline as it is 
progressively running out of gas. At the point where the 
cooling time overtakes the heating time, the two are of the 
order of ^ 1 Gyr. It is interesting to observe that although 
gas accretion onto the central galaxy has already stopped at 
z ^ 4, when Mhaio lO^^/i^^M©, the accretion of cold gas 
in substructures along the filaments continue to be impor- 
tant until z ^ 2, when Mhaio ~ 10^^h~^ Mq. Semianalytic 
models of galaxy formation explain the galaxy bimodality 
with the shutdown star formation in haloes above a criti- 
cal mass Merit- Our work partly explains why these studies 
find the best fit to t he data for Merit 2 x 10^^ M© (e.g. 
iGattaneo et al.ll2006l ) despite the fact that a purely theoret- 
ical analytic approach gives a shutdown mass for the disrup- 
tion of cold filaments of 3 x 10^^ M© (|Dekel BirnboimI 
I2OO6I V In conclusion. Fig. [6] and Fig. [7] show that, with- 
out refuelling by gas-rich mergers or cooling flows, the cen- 
tral galaxy is already on its way to getting 'red and dead' 
at z ^ 2.5. That occurs because shock heating rearranges 
the gas in the cold filaments into a hot static spherical 
halo (Fig. 2]). Without merging or cooling flows our galaxy 
would simply evolve into a red early-type spiral galaxy with 
Mstar - 1 - 2 X IO^^H-^Mq. 

Mergers cause galaxies to evolve toward earlier type 
morphologies but at the same time they reactivate star for- 
mation. In fact, mergers are the dominant mechanism in 
supplying cold gas to the galaxy after the accretion of cold 
filamentary flows has ceased. We caution that this may be 
partly due to our artificially long cooling time, which derives 
from computing radiative cooling for primordial metal abun- 
dances. Our cosmolo gical simulations supp ort the earlier 
conclusion derived bv lSpringel et al.l (|2005bl ) in the case of 
isolated mergers: mergers cannot form red elliptical galaxies 
unless there is a quenching mechanism that prevents them 
from leaving tails of star formation inconsistent with the 
red colours of elliptical galaxies (e.g. Springel et al. 20053)- 
Our simulations show that the model without AGN feedback 
fails to form a convincing elliptical galaxy. In the simulation 
without black hole, the galaxy ends up with a visual appear- 
ance and a Sersic index (n ^ 2.5) consistent with those of 
an elliptical galaxy. It also has boxy isophotes. However, it 
has a central blue light excess that is abnormal for a galaxy 
with these properties and that places it squarely on the blue 
sequence fFig. [TO)) . 

In contrast, already ^ 1% of the energy released by the 
growth of a supermassive black hole is sufficient to quench 
star formation in galaxy mergers if it can be coupled effec- 
tively with the surrounding gas. This claim is motivated by 
considering that our black hole has accreted a mass that is 
^ 5 times lower than the typical mass of the central black 
hole in an elliptical galaxy with Mstar 10^^h~^ Mq. There- 
fore, the same energy that we extract with p = 0.05 could 
have been extracted with a lower heating efficiency but a 
larger mass growth. The limited growth of the black hole is 
due to the promptness of AGN feedback, which immediately 
blows away the fuel for accretion as soon as the accretion 
starts. We believe that the lower growth in our cosmologi- 
cal simulation with respect to the growth of black holes in 
simulations of isolated mergers is partly a resolution effect, 
which causes the energy injection to be smoothed over a 
larger part of the galaxy. In simulations of isolated merg- 
ers, outflows are primarily affecting the central region, and 
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therefore they do not have the same immediate effect on gas 
at larger radii, which may be about to fall onto the black 
holes. 

Our simulation with AGN feedback forms a galaxy that 
is already red after z ^ 0.6 and that is perfectly on the red 
sequence at 2; = (Fig. llOp consistently with its elliptical 
morphology. The suppression of star formation, which has 
moved the galaxy to the red sequence, is due to the combina- 
tion of two physical phenomena: the shutdown of cold mode 
accretion by shock heating, and AGN feedback quenching 
star formation in galaxy mergers. AGN feedback also con- 
tributes to raising the entropy of the hot IGM by removing 
low entropy tails vulnerable to developing cooling flows. 

AGN activity appears to be directly connected with 
merging. The two times when the black hole reaches the 
Eddington limit {Mb ^ -22.6 and Mb ^ -23) coincide 
with the two major mergers at z ^ 2.04 and z ^ 0.075. 
We also find AGN activity without us being able to ob- 
serve an ongoing merger {z ~ 0.63), but, even in these 
cases, we see morphological disturbances and tidal features, 
which can be traced to the interaction with small compan- 
ions. It is in fact an episode of this kind that shuts down 
the residual star formation activity at z 0^ 0.63 and moves 
the galaxy from the blue to the red sequence. This episode 
supports the observational notion that AGN host galaxies 
are transition objects in migration from the blue to the red 
sequence (Graves et al. 2007; Schawinski et al. 2007), al- 
though we also see a transition in the properties of AGN 
hosts from blue and star-forming at 2; ~ 2 to red and dead 
at ^ - (Fig. [7|. This trend is in agreement with obser- 
vational results, which find that the hosts of low redshift 
AGNs are ma inly on the red sequence or on the top o f 
the blue cloud (iKauffmann et al.|[2003bl:ljahnke et allliooi a: 
ISanchez et "alll2004l : iNandra et al.ll2007l ). while at high red- 
shift the observations supp ort a connection b etween AGN 
activity and star formation (|jahnke et al ] |2004b ). The black 
hole accretion event at 2; ~ 0.1 is significant because it shows 
an example of an Mb ^ -23 AGN in a red (U - V ^ 1.2) 
boxy elliptical with a Sersic index of n :^ 6.8. One should 
not be puzzled by the fact that one of the two most impor- 
tant episodes of luminous accretion is happening at such a 
low redshift because we have simulated the formation of an 
intermediate mass elliptical in a halo that has just passed 
the critical mass for star formation shutdown. Therefore it 
is not surprising that this galaxy is hosting Seyfert activity 
at low redshift. Had we simulated the central galaxy of a 
galaxy cluster, where the critical halo mass is crossed much 
earlier, we would have seen a much more luminous AGN 
and a shutdown of AGN activity at a much higher redshift. 
This is directly related to the downsizing of galaxy formation 
(Cattaneo et al., in preparation). 

Our simulations contain metal enrichment of the ISM by 
supernovae. We have found that even moderate AGN winds 
can easily produce the required enrichment of the IGM at 
z ^ 3 although in realist ic situations supernova win ds are 
also likely to play a role (ISpringel Hernauistlliooi ). The 
field of studying the complex interaction between galaxy for- 
mation, AGNs and the ISM is still in its infancy. However, 
the need for AGN feedback in galaxy formation appears as 
a fundamental physical results independently of the many 
uncertainties concerning the details of how it works. 
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